Generalized similarity in finite range solar wind magnetohydrodynamic turbulence 
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Extended or generalized similarity is a ubiquitous but not well understood feature of turbulence 
that is realized over a finite range of scales. ULYSSES spacecraft solar polar passes at solar mini- 
mum provide in situ observations of evolving anisotropic magnetohydrodynamic turbulence in the 
solar wind under ideal conditions of fast quiet flow. We find a single generalized scaling function 
characterises this finite range turbulence and is insensitive to plasma conditions. The recent unusu- 
ally inactive solar minimum -with turbulent fluctuations down by a factor of ~ 2 in power- provides 
a test of this invariance. 
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A characteristic feature of fully developed turbulent 
flows is an inertial range which is scale invariant in its 
statistics. Classical signatures of this scale invariance 
in isotropic homogeneous hydrodynamic turbulence are 
power law power spectra with a Kolmogorov exponent 
of —5/3 and multiscaling in the structure functions of 
the velocity field. This scale invariance is modified when 
a characteristic scale becomes apparent, as is the case 
for measured flows at large but finite Reynolds number 
which are realized over a finite range of scales. There is 
longstanding interest in the nature of this correction to 
ideal hydrodynamic (Kolmogorov) turbulence for exam- 
ple in forced confined flows [l| and wall bounded shear 
flow [2| . The nature of this correction captures features 
of the structure and dynamics of the turbulent flow 0, Q 
and when seen from the perspective of finite sized scaling 
makes contact with critical phenomena 0] . A ubiquitous, 
but not well understood, aspect of finite range turbulence 
is generalised scale invariance or Extended Self Similarity 
(ESS) 6] which is seen in both hydrodynamic and magne- 
tohydrodynamic (MHD) turbulent flows, for example in 
the solar wind [7|, |_S|, |9j, [lfj . In this Letter we find that for 
solar wind MHD turbulence a single generalized scaling 
function captures statistical scaling of the largest (outer) 
scales and is at the heart of the observed ESS, and is 
insensitive to conditions of the plasma flow. This points 
to a universal property of finite range MHD turbulence. 

The solar wind provides a laboratory to study MHD 
turbulent fluctuations [Tlj in a high Reynolds number 
(l2| magnetofluid. In situ spacecraft provide single-point 
plasma observations over a range of scales from seconds 
to years. Typically, power spectra of solar wind velocity 
and magnetic field components show a power law scaling, 
suggesting the presence of an inertial r ang e of turbulence 
and at lower frequencies a "~ 1/ /" [l3l ITij spectrum. 
At higher frequencies there is a cross-over to a dissipa- 
tive or disp ersive range which may support a cascade 
[HI d, E3 EM 13- At greater solar distances, and in 
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slow as compared to fast solar wind, this inertial range 
extends to lower frequencies [201 ] consistent with an ac- 
tively evolving turbulent cascade. The inertial range tur- 
bulence is intermittent with non-Gaussian statistics at 
small scales [U [H [H, [H| and anisotropic [M IM E3 
due to a preferred direction introduced by the presence 
of a background magnetic field, and may also be com- 
pressible [28|,[29j]. Ensembles of fluctuations in the "1//" 
region are anisotropic [3(| but it is suggested that should 
conditions of isotropy hold locally, there is evidence for 
a turbulent cascade on these large scales [29|, |3l| . 

Observations in quiet, relatively uniform background 
magnetic field and flow over extended intervals are pro- 
vided by the solar polar passes of the ULYSSES space- 
craft at times of minimum solar activity. The most re- 
cent solar minimum has been reported to be anomalously 
quiet. Compared to the previous minimum, ULYSSES 
observations show a ~ —17% lower plasma density and a 
~ 15% lower magnetic field on average; the fluctuations 
in magnetic field are ~ 42% lower in power [H, [H, [HJ. 
We will compare ULYSSES polar passes for these two su- 
cessive solar minima to test the robustness of generalized 
scale invariance over different conditions. We compare 
10 day intervals of ~ 13,000 points each of one minute 
averaged magnetic field component measurements Bi(t) 
for each of four polar passes, the North and South passes 
for the two solar minima. The length of these data inter- 
vals is sufficient to explore the inertial range scaling and 
its transition to "1//" at lower frequencies. Contiguous 
intervals are selected when ULYSSES is deep within the 
polar coronal hole flow such that the heliospheric dis- 
tance is within [1.92 — 2.83 AU] and latitudes are above 
70° (and peak at 82°) [35|. The selected observations are 
of fast quiet polar flows with flow velocity ~ 750 km/s 
and are free from large transient coronal events. We have 
analyzed in total 24 intervals over the 4 polar passes and 
find that these all yield the same scaling behaviour in 
the inertial range that we report here. All results are 
in RTN coordinates: R is along the Sun-spacecraft axis, 
T is the cross product of R with the Sun's rotation axis 
and N completes the right-handed system; the ambient 
field is almost radial. Figure 1 compares a 10 day inter- 
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FIG. 1: Log-log plot of the R magnetic field component PSD 
for 10 day intervals from ULYSSES North polar passes in 1995 
and 2008. ULYSSES was at a similar heliospheric distance for 
these intervals. There is a difference in power of ~ —42% be- 
tween the two minima. Inset: Compensated GSF Sz/t (linear 
scale) versus r (logarithmic scale), for the R component for 
all passes at this heliospheric distance. 



val from each of the North polar passes for the sucessive 
minima at ~ 2.2 AU. The power spectral density (PSD) 
for both intervals shows a ~ —5/3 exponent power law 
region with a cross-over to ~ 1// at lower frequencies, 
the traces track each other but the overall power level can 
be seen to differ by a factor of ~ 2. We focus on gener- 
alized structure functions (GSF) of the field fluctuations 
S p =< \Bi(t + T) — Bi{t)\ p > where < ... > is an ensemble 
average over t and where in these conditions of fast flow 
we invoke the Taylor hypothesis and use time interval r 
as our measure of scale. For finite datasets with non- 
Gaussian Probability Density Functions (PDFs), these 
statistics can be affected by the presence of large outliers 
in the PDF tails; we have used the method in [3& HtJ to 
verify the robustness of our results against any bias in- 
troduced by this source of uncertainty. Inset in Figure 1 
we plot the compensated third order structure functions 
Sz/t. These curves have different amplitudes due to dif- 
ferent power levels in the fluctuations but intriguingly all 
roughly have the same functional form. 

Let us now examine the GSF in more detail. In the 
main panel of Figure 2 we plot S3 versus r for six consec- 
utive 10 day intervals. The inertial range extends from 
below the smalle st ( 1 minute) timescale up to r ~ 20 
minutes (see also [T3, H3| ) , and we can see that over this 
range there is no single straight line fit, as would be the 
case for fully developed turbulence for which one antici- 
pates the scaling S p ~ . Extended Self Similarity @ 
instead suggests the scaling of ratios of the GSF: 



FIG. 2: GSFs of S3 versus r on log-log axes for the R compo- 
nent of the magnetic field for 5 consecutive 10 day intervals 
from the North polar pass of 2008 shown for r = 2 — 60 min- 
utes. The inset figure shows ESS plots of S3 against Sa on 
logarithmic axes, the different intervals are shifted for clarity. 



and we test this in the inset of Figure 2 where we plot 
S3 versus S2. We can see that ESS indeed captures the 
scaling of the inertial range (but does not extend into the 
"1//" range, for a detailed ESS analysis of the 1995 polar 
pass see QJl). We find that ESS holds for the inertial 
range of all of the intervals under study. ESS implies a 
generalised similarity of the form: 



S p (r)=S p (T )[g(r/r )}^ 



(2) 



however the observation of ESS alone does not discern 
whether or not the function g(r/ro) is ubiquitous in na- 
ture or whether it varies across different components of 
the magnetic field or with ambient conditions. We first 
establish that the same function g is obtained for differ- 
ent orders p (as required by ESS) and also for all three 
components of the field and for each of the polar passes. 
To obtain an expression for <?(t/to) we need to invert 
equation (2) however since we do not have the scaling 
S p ~ t^P- 1 we cannot determine the exponents di- 
rectly. Instead we rearrange equation (2) to give: 



9(r,T ) 



S P (r) 



S p (t ) 



OA) 



= g(r/ro) Ci3) 



(3) 



S p (t) ~ [S q (r)} 



C(P)/C(9) 



(1) 



ESS plots of one GSF against another (as in the inset 
of Figure 2) then yield measurements of the ratios of 
the exponents ((3)/((p). We plot g versus r on linear 
axes in Figure 3 for orders p = 1 — 4 and for all of the 
magnetic field components. Each of the four panels refers 
to a 10 day interval from each of the polar passes, again 
at the same heliospheric distance. The plotted range of 
t = 1—60 minutes encompasses the inertial range and the 



FIG. 3: [S p (t)/S p (t )] c{3)/c{p) versus t/to, (to = 10 minutes) 
for p = 1 to 4 for all field components are overplotted for 10 
day intervals at the same heliocentric distance for each of the 
ULYSSES polar passes: South 1994 (days 260 - 269), North 
1995 (days 240 - 249), South 2007 (days 60 - 69) and North 
2008 (days 40 - 49). 



cross-over to the "1//" range. The normalization scale To 
is arbitrary and we used To = 10 minutes which is within 
the inertial range for all the intervals. In all four panels 
of Figure 3 we can see two distinct regions: an inertial 
range at smaller t where the curves closely correspond 
to each other, and the "1//" range at larger r, where 
the curves diverge. Variability in the scaling exponent a 
in the T// Q " , a ~ 1 region has been noted previously, 
both with ambient conditions [3(| [H, Hil and with field 
component Figure 3 demonstrates that there is a 

robust (i.e. usually present) inertial range which follows 
a generalized scaling captured by function g(r) which 
satisfies ESS (i.e. it is the same for allp). Intriguingly, it 
is also the same for all of the field components, despite the 
differences in their power levels, and we have verified that 
this is the case for all of the intervals under study. The 
outer scale t c where the inertial range terminates varies 
between all four polar passes and is not simply ordered by 
the distinct ambient conditions of the two solar minima; 
the longest (t c ~ 25 minutes) and the shortest (t c ~ 12 
minutes) both occurring respectively for the 1994 and 
2007 South solar minima. 

We next test whether the same function g(r) holds for 
all of the polar passes. Figure 4 again plots g(r) versus r 
on linear axes, now for p = 3 for all field components for 
all four polar passes shown in Figure 3. Again we find an 
inertial range where there is a close correspondence be- 
tween all these curves. A good fit to this function is of the 
form h(r,To,f) — a.(r/f )( T / r °) where we have chosen 
scale parameters To = 10 minutes and f = 1 minute and 
obtained constant fitted parameters a = 0.101±0.001 and 
b = 0.10 ±0.01, this is overlaid on the figure. This single 
function g reflects both the underlying phenomenology 



FIG. 4: [S 3 (t)/S , 3 (to)] c(3)/c(3) versus t/t for all field com- 
ponents for the 10 day intervals plotted in Figure 3. 



of the inertial range MHD turbulence in the value of the 
exponent £(3) and the (generalized) scale invariance of 
the largest structures in the evolving turbulence in the 
function g. To recover the same g(r) over this variety 
of ambient conditions suggests invariant phenomenology 
for both these aspects of this finite range, evolving MHD 
turbulence. Although we cannot precisely determine £(3) 
the robustness of the inertial range would suggest that it 
does not change significantly between the different inter- 
vals. If so, we have shown an important property of MHD 
turbulence as seen evolving in the solar wind- that ESS 
arises as a consequence of a single generalized similar- 
ity of finite range turbulence that is invariant to ambient 
conditions and the level of power in the turbulent signal. 

To gain some insight into the robustness of this gener- 
alized similarity we finally turn to the PDFs of the fluctu- 
ations Bi (t +t) — Bi (t) . Figure 5 shows the PDFs for the 
radial field component fluctuations in the inertial range 
(t = 6 minutes shown) for all of the intervals in Figure 
3, normalized to a to afford a comparison of the func- 
tional form. We can see that all of the curves collapse 
onto each other - the functional form of the turbulent 
fluctuations is invariant to different power levels in the 
turbulence and different ambient conditions. The PDF 
is non- Gaussian (a Gaussian is fitted for comparison), a 
well-known feature of small scale solar wind turbulence, 
and is stretched exponential in form. All three compo- 
nents show this invariance to ambient conditions in the 
functional forms of their PDFs. The field components 
differ in the power levels of their fluctuations and this is 
captured by the standard deviation a. For the R com- 
ponent cr = [0.37 - 0.3]nT ± 0.002 whereas for the T and 
N components a = [0.52 - 0.42] raT ± 0.003 between the 
1994/5 and 2007/8 solar minimum passes. For compar- 
ison, the mean magnetic field is ~ 1.4 nT (94/95 polar 
pass) and ~ 1 nT (07/08 polar pass). There is therefore a 
consistent anisotropy in the amplitude of fluctuations be- 
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FIG. 5: The PDF of fluctuations in the R magnetic field 
component on timescale r = 6 minutes for 10 day intervals 
plotted in Figure 3. A Gaussian fit to the data is also plotted. 

tween the R and the T, N components, present for both 
solar minima as well as an overall decrease in power be- 
tween minima. 

The recent unusually quiet solar minimum provides a 
unique opportunity to make a relatively controlled com- 



parison of the statistical scaling of evolving MHD turbu- 
lence in situ in the quiet, fast solar wind under different 
ambient plasma parameters (density, magnetic field mag- 
nitude) and in the power level of the turbulence. We find 
evidence of the same generalized scaling in all 24 inter- 
vals that we consider. The robustness of this generalized 
similarity suggests that this may be a universal property 
of finite range anisotropic MHD turbulence. How this 
relates to the statistical scaling properties of the evolv- 
ing solar wind under all manifestations- such as in slow 
as opposed to fast flow, and over the solar cycle (37j - 
is an open question. Opportunities to test for this gen- 
eralized similarity in finite range turbulence are also in 
the turbulent foreshock [Z(| and in the magnetosheath 
[4ll ] provided sufficiently long stationary intervals can be 
identified. 
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